Surface Critical Behavior of bcc Binary Alloys 
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The surface critical behavior of bcc binary alloys undergo- 
ing a continuous B2-A2 order-disorder transition in the bulk is 
investigated in the mean-field (MF) approximation, employ- 
ing a semi-infinite lattice model equivalent to an Ising anti- 
ferromagnet in an external field. Our main aim is to present 
clear evidence for the fact that surfaces which break the two- 
sublattice symmetry generically display the critical behavior 
of the normal transition, whereas symmetry-preserving sur- 
faces exhibit the behavior of the ordinary transition. To this 
end, the lattice MF equations for both symmetry-breaking 
(100) and symmetry-preserving (110) surfaces are cast in the 
form of nonlinear symplectic maps, the associated Hamilto- 
nian flows are analyzed, and the length scales involved are 
computed. Careful examination of the continuum limit yields 
the appropriate semi-infinite Ginzburg-Landau model for the 
(100) surface and reveals subtleties overlooked in previous 
work. The continuum model involves an "effective" order- 
ing surface field g\ 7^ 0, which depends on the parameters 
of the lattice model. The singular behavior predicted by the 
Ginzburg-Landau model is shown to agree quantitatively with 
the solutions of the lattice MF equations. 

68.35.Rh, 64.60.Cn, 05.50.-Hq 



I. INTRODUCTION 

Experiments on binary {AB) alloys that undergo an 
order-disorder transition in the bulk have yielded a 
wealth of information on surface critical phenomena in 
semi-infinite matter.^ In these systems one inevitably has 
to cope with the influence of surface segregation, i. e., 
the enrichment of one component at the surface. Sur- 
face segregation occurs, e. g., due to different interaction 
energies or sizes of the two species. Theoretically, the 
variation of the local composition near the surface may 
necessitate the introduction of "nonordering" densities, 
which are given by linear combinations of the local con- 
centrations of A and B atoms on the various sublattices. 
In the case of surface critical phenomena at first- order 
bulk transitions, such as "surface-induced disordering" 
in fee binary alloys, nonordering densities strongly influ- 
ence the asymptotic behavior.^ In this paper we are con- 
cerned with bcc alloys that exhibit a continuous (second- 
order) bulk transition and are thus promising candidates 
for testing current theories of surface critical behavior at 
bulk critical points. '^"^ 

The continuous B2-A2 transition occurring in FeAl or 
FeCo has been investigated previously by F. Schmid.^ 



She studied a semi-infinite lattice model equivalent to 
a bcc Ising antiferromagnet both by Monte Carlo sim- 
ulation and within the mean-field (MF) approximation, 
and made the important observation that the orienta- 
tion of the surface in general matters. Her conclusions 
can be summarized as follows: (a) A nonvanishing order 
parameter (OP) profile occurs for T > Tc, the bulk criti- 
cal temperature, provided that (i) the surface breaks the 
two-sublattice symmetry (see below), and (ii) one com- 
ponent is enriched at the surface, (b') The observable 
surface critical behavior should be representative of the 
ordinary universality class even if the above conditions 
(i) and (ii) are met. 

While we agree with (a), we find that (b') should be 
replaced by (b): If conditions (i) and (ii) are satisfied, 
the surface critical behavior generically is characteristic 
of the normal transition, which belongs to the same uni- 
versality class as the extraordinary transition.^ 

In a foregoing Letter* by Drewitz, Burkhardt, and our- 
selves, exact transfer matrix calculations were employed 
in conjunction with conformal invariance to present clear 
evidence for (a) and (b) in bulk dimension d = 2. Here we 



(a) 



(110) surface 




(b) 



(100) surface 




FIG. 1. Examples of (a) symmetry-preserving, and (b) sym- 
metry-breaktng surfaces. Sites of sublattices a and /3 corre- 
spond to open and full circles, respectively. In the A2 phase, 
the concentration of either component is the same on all sites, 
whereas the two sublattices are preferentially occupied by A 
and B atoms, respectively, in the ordered B2 phase. 
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generalize these results to arbitrary d using MF theory 
and a mapping onto a Ginzburg-Landau model. 

The reason for the appearance of normal critical be- 
havior is a subtle interplay between the symmetry with 
respect to sublattice ordering and broken translational 
invariance due to the free surface. Consider first a finite 
system with periodic boundary conditions. The precise 
form of the Hamiltonian Ti.{(7i} does not matter here 
and will be given in Sec. 11. The spin variable cji — 1 
((Ti = — 1) represents an A (B) atom on lattice site i. 
The statistical weight of a configuration {ai} is given by 
the finite- volume Gibbs measure 

piWi}) ^ ^e-^^'^'^'K P^l/{kBT), (1) 

where T is the temperature and ks denotes Boltzmann's 

constant. The normalization factor Z is the grand- 
canonical partition function. The local concentration Ci 
of A atoms at site i can be expressed in terms of the mean 
magnetization {(Ti) as Ci = (1 + {ai))/2. The Gibbs mea- 
sure (1) is translationally invariant, 

PiWi}) = PiWi}) , '^i = <^i+t , (2) 

where t may be chosen arbitrarily from the set T of 
all bcc lattice vectors. Due to spontaneous symmetry 
breaking, the thermodynamic states that are obtained by 
calculating expectation values with the measure (1) and 
taking the infinite- volume limit need not be translation- 
ally invariant. (Formally, one introduces a symmetry- 
breaking "staggered" field which is sent to zero after the 
thermodynamic limit has been performed.) In the or- 
dered phase below Tg, one has {ai) ^ {<^i+t) for all trans- 
lations t = taf3 that map a to /3-sites (Fig. 1). 

Surfaces are introduced by imposing free boundary 
conditions along one direction while retaining periodic 
boundary conditions in the other directions. Then the 
measure (1) is invariant under the subset T' C T of 
translations parallel to the surface. We call the sur- 
face symmetry-preserving if T' contains a "sublattice- 
exchanging" translation t^^, and symmetry-breaking oth- 
erwise. In the case of the bcc lattice considered here, 
symmetry-breaking surfaces arc characterized by the al- 
ternation of a and /3-planes along the direction normal 
to the surface. Let us assume that no spontaneous sym- 
metry breaking takes place above Tc, which would re- 
quire supercritically enhanced surface couplings, so that 
{oi) = {(Ti+t) for T > Tc and t e T'. 

Thus for symmetry-preserving surfaces and T > Tc, 
both sublattice magnetizations are the same within each 
plane parallel to the surface, and the OP profile vanishes. 
Nonetheless one obtains an inhomogeneous magnetiza- 
tion (or concentration) profile due to surface segregation. 
E. g., segregation of one component at the symmetry- 
preserving (110) surface leads to an alternation of A and 
B-rich lattice planes since the interactions favor the occu- 
pation of nearest-neighbor (NN) sites by different species. 
As will be shown in Sec. Ill A, the scale on which the seg- 
regation profile decays to its bulk value is on the order of 



the lattice constant even near Tc- This length should be 
irrelevant (in the renormalization-group sense), so that 
the surface must exhibit ordinary critical behavior. 

Symmetry-breaking surfaces, like the (100) surface, de- 
stroy the two-sublattice symmetry. Surface segregation 
again leads to an inhomogeneous concentration profile for 
T > Tc, which is now equivalent to a nonvanishing OP 
profile since adjacent lattice planes belong to different 
sublattices. The OP profile decays on the scale of the 
bulk correlation length, which diverges for T ^ Tc (Sec. 
Ill A). Such a persistence of surface order for T >Tc has 
been confirmed in recent experiments on FeCo(lOO).® 

According to the experimental results of Ref. 9, su- 
percritical enhancement of the surface couplings can be 
ruled out. Thus it is natural to attribute the persis- 
tence of surface order to an ordering surface field gi ^ 0. 
This field causes the system to display the normal tran- 
sition. However, for the binary alloys considered here an 
ordering field corresponds to a local chemical potential 
acting differently on the two sublattices (staggered field 
in magnetic language). There is no natural source for 
such a field on the microscopic level. The challenge is 
to demonstrate in an unequivocal fashion that a nonzero 
gi nevertheless emerges in a continuum (coarse-grained) 
description, i. e., in the context of a Ginzburg-Landau 
model, and to derive a MF expression for gi in terms of 
the lattice model parameters. Of course, in comparing 
theory and experimental or simulation data, one should 
keep in mind that gi may be small, so that the crossover 
to normal critical behavior occurs only close to T^.^^ 

In the next section, we shall reformulate the lattice MF 
equations for the semi-infinite alloy with free (100) and 
(110) surfaces as a problem in discrete dynamics, i. e., 
the iteration of nonlinear symplectic maps.'^"'^ From the 
linearized maps the characteristic length scales of both 
the concentration and OP profiles will be calculated (Sec. 
Ill A). The full nonlinear maps will be analyzed in Sec. 
Ill B. After the derivation of the Ginzburg-Landau model 
for the (100) surface in Sec. IV, we shall compare the 
predictions of the continuum theory with the numerical 
solutions of the lattice MF equations (Sec. IV B). This 
will serve to demonstrate that the (100) surface displays 
normal critical behavior generically. However, for cer- 
tain exceptional parameter values, where gi happens to 
vanish at T = Tc, one recovers the singularities of the 
ordinary transition (Sec. IV C). We will summarize our 
main results in Sec. V. 



II. MF EQUATIONS AS NONLINEAR MAPS 

We consider the lattice-gas model of a binary (AB) 
alloy on a bcc lattice. Each atomic configuration is char- 
acterized by the values of the occupation variables pf, 
pf , where p'- = 1 if site i is occupied by an atom of type 

E {A, B} and p^ = otherwise. Within the grand- 
canonical ensemble, the configurational energy reads 
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i^j V,T V i 

where e^J = eJJ is the interaction energy between v and 
T atoms at sites i and j, and /i^i and /xs are chemical 
potentials for A and i? atoms, respectively. We neglect 
vacancies, so that pf +Pi = 1, and rewrite the occupa- 
tion variables in terms of Ising spins = ±1 as 

P^ = ^(l + ^i), pf = \{l-ai). 
Then (3) takes the form of an Ising Hamiltonian 

^Wi] = E '^ij'^i^i - E ^^^^ ' '^^^ 

tT^i i 

where a spin-independent term has been dropped and 
1 



= i -eft- ' (4a) 

= I E (^S"" - 4^) + J (/"A - Ms) • (4b) 
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In the following we only consider NN interactions e"^"*, 
e^^, and e^^. Moreover, we do not allow for enhanced 
surface interactions. For an B2-A2 order-disorder transi- 
tion to exist , the Ising coupling J = (2e'^^ - e^^ - e^-^ ) /4 
must be antiferromagnetic (J < 0). For semi-infinite 
systems with (100) or (110) surfaces the local field (4b) 
differs from its bulk value only in the first layer, 



Hi 



where 



H + Hi if i e surface , 
H otherwise , 

if=^(e^^-e^^) + i(MA-M, (5a) 



H, = 



Cs-C 



(^ 



BB AA 



(5b) 



Here, ( and (^s are the coordination numbers of bulk and 
surface spins, i. e., C = 8, while C« = 4 and Cs = 6 for the 
(100) and (110) surface, respectively. 

Some remarks about the role of the surface field Hi 
are in order here. The field Hi favors one component at 
the surface and thus accounts for surface segregation (see 
Sec. I). Because of different interaction energies e^"^ ^ 
e^^ it is nonzero gencrically. More generally. Hi also 
models other effects such as different sizes of the two 
constituents. For symmetry-preserving orientations. Hi 
acts uniformly on a and /3-sites at the surface and must 
not be confused with an ordering (staggered) field. For 
symmetry-breaking surfaces, spins on a and /?-sites in 
the first two layers experience different fields H + Hi 
and H. Hence Hi should contribute to an "effective" 
staggered surface field gi ^ 0. However, even if Hi = 



(but H ^ 0) one obtains an inhomogeneous (oscillating) 
magnetization profile equivalent to a nonzero local order 
parameter, and an ordering surface field gi ^ should 
again emerge in a coarse-grained description. 

The mean-field (MF) or Bragg- Williams approxima- 
tion is conveniently formulated in terms of a variational 
principle. The free-energy functional reads 



-FMFA{(o-i)} = W{(cri)} - TSuFAiicTi)} , 

with the MF entropy 

'S'MFA{(ci)} = —ks 2. / d,^ tanh x . 
i Jo 

Variation of FMFA{(fi)} yields the MF equations 



(6) 



((Ti) = tanh 



[Hi + ^Jijiaj) 



kBT 



(7) 



The magnetization densities of the two siiblatticcs vary 
only in the direction perpendicular to the surface. For 
the (100) surface one may thus write 

((Ti) = m„ for i G lattice plane n , 

where m„ is the magnetization density of lattice plane n. 
Likewise one has, for the (110) orientation. 



mn,a for i € plane n, sublattice a , 
rrin^is for i G plane n, sublattice P . 



(5) It is convenient to introduce the reduced quantities 



K 



Then the MF equations read, for the (100) surface, 

^ tanh"^ mn = h- m„_i - m„+i 
K 



for n > 1, and 

tanh~^ mi = h + hi — 1712 ■ 
For the (110) orientation, one has 



— tanh mn,a = h — m. 



(8) 



(9) 



(10a) 



— tanh mn,i3 = h- m„,c< ■ — ^ — (10b) 



for n > 1, and 



--- tanh ^ TOi a = h + hi 



K 



tanh ^ fni,i3 = h + hi — mi^a — ^fn2,a ■ 



(11a) 
(11a) 
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Wc now combine the magnetization densities of two 
neighboring planes into single points in and M^, re- 
spectively, 

v„ = (m„_i,m„)^ , 

Wn = (m„_i,a,m„_i,/3,m„,a,'Tln,/3)'^ , 

where "T" denotes the transpose. Then (8) and (10) are 
equivalent to the recursion equations 

Vn+l = F{Vn) , (12) 
Wn+1 = G{Wn) , (13) 

where the nonlinear maps F and G are defined by 



where ma, rnf} are the solutions of the bulk MF equations 
[Eqs. (A2)]. As will be discussed in the next section, the 
bulk solutions correspond to fixed points of the maps F 
and G. Then (18) implies that the trajectories converge 
to these fixed points. The task of solving the MF equa- 
tions for the semi-infinite system (for given values of K, 
h, and hi) thus translates into finding the intersections 
of the stable manifold (inset) of the corresponding fixed 
point with the surface boundary conditions (16) and (17) 
(see Sec. IIIB). 

We finally quote an important symmetry property of 
the above maps. The MF equations (8) and (10) are 
symmetric with respect to interchanging the layer mag- 
netizations of the planes n — 1 and n+1. Thus one has 



and 



G 



F : 







: 





tanh ^ ; 



V 



X3 

2h- 
2h- 



\ 



Xi — 2x3 — ^ tanh ^ X4 
X2 — 2x4 — ^ tanh~^ X3 y 



One advantage of rewriting the MF equations in terms of 

the discrete dynamics (12) and (13) is that one may gain 
an overview of all solutions by iterating arbitrary start- 
ing points Vi and Wi. In this way one obtains trajecto- 
ries Vi,V2, ■ ■ ■ and Wi, W2, ... in a two-dimensional (2D) 
and four-dimensional (4D) phase space, respectively. The 
maps F and G are both symplectic, i. e., their differentials 
DF and DG are symplectic matrices, and thus generate 
a discrete Hamiltonian dynamics on these phase spaces. 
Note that any symplectic map is volume-preserving, in 
particular. The theory of nonlinear dynamics offers con- 
venient tools to understand the discrete dynamics gener- 
ated by such maps."'^'' 

The MF equations (9) and (11) take the same form 
as (8) and (10) if fictitious zeroth layer magnetizations 
mo = —hi and mo, a = Jtiq,/? = —2hi are introduced, i. e., 

v^ = F{vi), (14) 
w:, = G{wi), (15) 

where Vi and Wi satisfy the boundary conditions 

vi = {-hi,mif, (16) 
wi = {-2hi, -2hi,mi^oi,rni^pf ■ (17) 

Moreover we require that the sublattice magnetization 
densities approach their bulk values for n — > 00, 



m„ 



m2n 



1710. 



(18) 





( m„+i 














I rrin , 




i TO„_i j 



and an analogous relation for G. It follows that both 
maps are invertible and that their inverses and G~^ 
are given by 



(19) 



where 



i^-i = RoFoR, 
G-'^ = SoGoS, 



R : (a;,2/)^ {y,xf , 

S : {xi,X2,X3,X4Y' ^ {X3,X4„Xi,X2)'^ ■ 



III. ANALYSIS OF NONLINEAR MAPS 
A. Linearized MF equations and length scales 

The bulk MF equations (A2) arc equivalent to the fixed 
point equations of the nonlinear maps G and = FoF, 
the second iterate of F. For T ^Tc, linearization of the 
maps about the fixed points yields the asymptotic (ex- 
ponential) decay of the sublattice magnetization profiles 
away from the surface. The decay lengths can be ex- 
pressed by the eigenvalues of the linearized maps. In the 
case of the (100) surface only one length scale, propor- 
tional to the OP correlation length, governs the decay of 
both the sublattice magnetization and the OP profiles. 
An additional length associated with the mean magneti- 
zation density occurs for the (110) orientation. However, 
this length scale stays finite at T = T^. 

An analysis of the bulk MF equations is straightfor- 
ward and may be found in Appendix A. The critical 
coupling Kc = {4\J\) / {ksTc) as a function of the uni- 
form bulk field h is determined by 



1 



2Kc{h) 



= 1 - mdhf 



(20) 
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where rric = mricQi) is the uniform magnetization at T = 
Tc- Using (A6) to ehminate rric in favor of Kc and h in 
(20), one obtains an expression for the critical hne (Fig. 

2)14,15 

The bulk sublattice magnetization densities may be 
written 



m„ = m 



mp — m — (/) , 



(21) 



where m and (j) are the mean magnetization and the OP, 
respectively. In the region of the phase diagram where 
the A2 phase is thermodynamically stable, the bulk MF 
equations have a unique solution (see Appendix A) 

m^mdis{K,h), 0^=0. (22) 

Thus the only fixed points of and G in this case are 

\T 



VdiB = (mdis,"idis) 
Wdis = {rridis, "idis, mdis, "idis)"^ • 



(23) 



The solution (22) becomes thermodynamically unstable 
on crossing the critical line (Fig. 2). At the same time, 
two new solutions describing the pure B2 phases emerge, 

m^m,,diK,h), (b = ±MK:h)^0, (24) 

corresponding to the fixed points 

1,2 



'^ord = (™ord T (t'b, "T'ord ± </>&) , 
■"^ord = (™ord ± (t>b, "^OI■d T 'f'b, "^ord ± 06, "^OI■d T ' 



(25) 



Note that the fixed points t'^d' ''^ord form a 2-cycle of 
the map F, 



2.5 





A2 phase 






B2 phase \ 
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V, , 


J 


1 



^ 1.5 

0.5 


-3-2-10123 

h = H/m) 

FIG. 2. MF phase diagram of the NN Ising antiferromag- 
net on a bcc lattice, showing a line of continuous transitions 
between the disordered (A2) and the ordered (B2) phase. 
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FIG. 3. Temperature dependence of the bulk OP and mean 
magnetization density. Below Tc, the disordered state (dashed 
line) becomes thermodynamically unstable. 



F (^ord) = -^ord , 



Flv 



ord ) 



^ord ■ 



(26) 



The bulk OP and the mean magnetization density at 
fixed bulk field h — 1 are shown in Fig. 3 as a function 
of the reduced temperature 



t 



T-T, 



K 



The OP vanishes as for t 0-, with the usual MF 
exponent (3=1/2 [Eq. (A8)]. The mean magnetization 
density is a nonordering (or noncritical) density. Typi- 
cally, such quantities display singularities, where a 
is the bulk specific heat exponent. In MF theory this be- 
havior reduces to a discontinuity in the first temperature 
derivative since a = [see Eqs. (A7) and (A8)]. 

The linearization of about any of the fixed points 



Vdis and V', has the eigenvalues 



Al,2 = -1 + 2UaUi3 ± 2^UaU(3y/uaUp - 1 

and the associated eigenvectors 



1 + Ai,: 
2w/3 



(27) 



(28) 
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T >T, 



IniA 



ImAf 




T <T, 



ImA,: 



ImA{ 




, ord 











A ord 


^ord^^ord 


1 dis 
^3 


/ ReA, 

_--^A dis 

^^2 



ReA, 



FIG. 4. Behavior of the eigenvalues (27) and (29) in the 
complex plane as one crosses the critical line (cf. Fig. 2) . The 
superscripts "dis" and "ord" indicate the linearizations about 
the fixed points (23) and (25), respectively. 



where Ua and Ujj are defined by (A3). Since F is symplcc- 
tic and thus area-preserving (Sec. II), one has A1A2 = 1. 

Likewise, the eigenvalues of the linearization of G about 

1,2 

tUdu 



Al,2 = -1 + 2^UaUi3 ± 2^UaU[) ~ ^UaUf3 . 
A3,4 = -1 - 2^UaU(i ± 2^UaUj3 + ^UaU0 , 



(29) 



with the corresponding eigenvectors 

Lin = I 1 



) Ai 2,— ^/ Ai 2 

Uf3 V U/s 



-£^3,4 = ( 1, — ,A3.4, ,/ — A3,4 



(30) 



Since G is syniplcctic, the eigenvalues come in pairs (Ai, 
A2) and (A3, A4) with A1A2 — A3A4 = 1. Repeated ap- 
plication of the linearized map to Li,2 generates eigenso- 
lutions with opposite sublattice magnetization densities 
within each layer. Therefore, Li^2 represent "ordering" 
eigenmodes. The eigensolutions generated by ^3,4 show 
an oscillating profile of the magnetization density because 
A3,4 < 0. However, the local magnetization is the same 
for all sites in a given lattice plane parallel to the surface. 
Hence the OP profile vanishes and one may refer to £3,4 
as "nonordcring" eigenmodes. 

The behavior of the eigenvalues in the complex plane 
as one crosses the critical line is shown schematically in 



Fig. 4. The eigenvalues Ai, A2 collide at +1 and form 
a complex conjugate pair on the unit circle in the or- 
dered phase. Thus the character of the fixed point v^is 
changes from hyperbolic to elliptic. At the same time, 
two new real eigenvalues corresponding to the hyperbolic 
fixed points v^^^ emerge. Within the notions of nonlin- 
ear dynamics, the map F undergoes a period-doubling 
bifurcation. The eigenvalues Ai, A2 of the 4D map G 
show an analogous behavior. However, the fixed point 
tUdis remains unstable in the ordered phase since A3, A4 
stay real. 

The solutions of the linearized MF equations satisfying 
the bulk boundary conditions (18) read 

V2n+1 =V.^.+ 02X2^ , 

Wn+1 =w^ + A2K^L2 + AzK'lLs , 

where and to* stand for one of the fixed points (23) 
and (25). The coefficients 02, A2, and are fixed by 
the surface boimdary conditions (16) and (17). 
Thus for the (100) orientation we obtain 



m2n+i =ma+ ae 



-(2n+l)/€ 



(31) 



»Ti2n =m^-be ^"/^ , 
with the decay length 

C = 2|lnA2r\ (31a) 
and the amplitudes 

a = {hi + rriis).— , b = hi +1710. (31b) 

V 

In the disordered phase, the amplitudes simplify to 

a = b = hi + rridis • (31b') 
We define the local OP (?i„ by 



i(-l)"(m„+i 



(32) 



where the power of —1 ensures that one always subtracts 
the magnetization densities of /3-planes from those of 
a-planes.^^ Eqs. (31) imply a nonvanishing OP profile, 
which decays on the scale of ^. In fact, ^ may be iden- 
tified, up to a proportionality factor, with the bulk OP 
correlation length. To sec this, one expands 2iaUp in Eq. 
(27) in powers of t using (A7) and (A8). This gives 



-1/2 



+ 0{t) , 



where 

C+ = e+(/i) 

e 



2y'2 - 2Kc{h)mc{h)[h - 2mc{h)] ' 



(33) 



(33a) 



(33b) 
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The decay length displays a \t\~'^ singularity just as the 
bulk correlation length, with the MF exponent — 1/2. 
Asymptotically, ^ should thus be proportional to the cor- 
relation length. Indeed, one has ^+{h)/^-{h) = v^, 
which is the MF value of the universal amplitude ratio of 
the correlation lengths above and below Tc-^"^ As t — > 0, 
the exponential decay of the OP profile becomes a power 
law, whose precise form will be investigated in Sec. IV. 

Likewise, the magnetization profiles for the (110) ori- 
entation are 



= ma + yle-"/«' + (-l)"ie-"/« , 
,0 = mi3- Be-"/^' + (-l)"Be-"/«, 



where now two length scales appear, 
C' = |lnA2r\ |=|ln|A3| 
The amplitudes are given by 



(34) 



(34a) 




(34b) 



In the disordered phase above Tc, they simplify to 

A = B = 0, A = B = -2hi - mdis • (34b') 
In particular, the OP profile 

= ^{'mn,a - mn,^) (35) 

vanishes for T > which is a consequence of the sym- 
metry of the (110) surface with respect to the two sublat- 
tices and the fact that neither enhanced surface couplings 
nor a staggered surface field are present. Asymptotically, 
^' and ^ behave as 



In 



(3 - 



0(t), 



(36a) 
(36b) 



with ^± given in (33) and |ln (3 - 2y/2) \ ^ ~ 0.57. The 
length ^' associated with the ordering eigenmodes di- 
verges ast and may again be identified with the bulk 
correlation length (up to a proportionality constant), "'^^ 
whereas £, stays on the order of the lattice constant. As 
will be seen in the next section, ^ describes the decay of 
the mean magnetization profile for T > Tc- 

Note that for both surface orientations the layer mag- 
netizations oscillate about mdis for T > Tc due to the an- 
tiferromagnctic coupling between adjacent lattice planes. 
However, only in the case of the (100) orientation does 
this oscillating profile lead to a nonvanishing OP profile 



whose characteristic length scale diverges as i — > 0+. In 
view of the presumed absence of enhanced surface cou- 
plings, such a behavior should be due to an "effective" 
ordering surface field g\ = g\{K,h,h\). Away from Tc 
and in the disordered phase, such a field generates a lin- 
ear response of the local OP which decays exponentially 
into the bulk. A glance at (31) leads us to anticipate the 
form 

9i(K,h,hi) ^ hi + mdis {K, h) . 

We will derive a formula for gi identical with the above 
expression in Sec. IV, when we map the lattice model 
onto a continuum theory. There it will become clear 
that gi is indeed a surface field coupling to the local 
OP that enters into a coarse-grained (Ginzburg-Landau) 
free-energy functional. 



B. Solutions of the nonlinear recursive maps 

The thermodynamically stable solutions of the bulk 
MF equations correspond to hyperbolic fixed points of 
the maps and G (cf. Sec. Ill A). The stable (un- 
stable) manifold Wg (VK„), or inset (outset), of a hy- 
perbolic fixed point consists of all points that converge 
to the fixed point under iteration of the map (inverse 
map). In order to solve the MF equations for the 
semi-infinite system one must determine the intersec- 
tions of the stable manifolds with the linear subspaces 
defined by the surface boundary conditions (16) and 
(17), i. e., the fine {{—hi,x)^\xG'M.} and the hyper- 
plane {(— 2/ii, — 2/ii, X, y)-^| x,y G M}. The magnetiza- 
tion densities of the first layer can be read off from 
the intersection points Vi = (— /ii,mi)'^ and Wi = 
{—2hi, — 2/?,!, mi^ct, mi^^)-^. The complete magnetization 
profile follows from the trajectories passing through vi 
and Wi.^^ Below Tc, infinitely many intersections exist 
and one has to resort to the original variational principle 
to find the equilibrium profile minimizing the free-energy 
functional (6). 



1. (100) surface, T>Tc 

In the disordered phase, the only fixed point of F"^ 
(as weh as of F) is t>dis = (mdis,mdis)^ [Eq. (23)]. Fig. 
5 shows a plot of the invariant manifolds for particular 
values of K and h. Fig. 6 depicts the magnetization 
and OP profiles obtained from the intersection with the 
boundary condition for hi = 0.4. Since Ro F = F~^ o R 
[cf. Eq. (19)], the inset and outset are mapped onto each 
other under reflection at the line x = y. 

The picture of the invariant manifolds does not change 
qualitatively as one varies the parameters K and h. Prom 
the eigenvector I2 [Eq. (28)] one infers that the slope of 
the inset at f dis vanishes in the high-temperature limit 
{K ^ 0), whereas it approaches —1{oj:K'\Kc{TI Tc). 
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FIG. 5. Invariant manifolds (inset and outset) of the hy- 
perbolic fixed point t^dis- The direction of the fiow under 
iteration of the map F is indicated by arrows. The dashed 
line represents the boundary condition for hi — 0.4. 



The upshot is that we find a unique intersection of the 
boundary condition with the stable manifold. In partic- 
ular, we obtain a nonvanishing order parameter profile 
at any temperature T > Tc. This conforms with the idea 
that the OP profile is due to an ordering surface field. 

An exceptional case occurs if the boundary condition 
exactly hits the fixed point Vdis, so that m„ = mdis and 
(j)n = 0. However, except in the case h = hi = 0, this 
can only be achieved for a special temperature T = Tq 
(at fixed h and hi). For T ^ To the OP profile is still 
nonzero. 



2. (110) surface, T > Tc 

As we have seen in Sec. I, this type of surface is 
symmetry-preserving and the Hamiltonian (4) is exactly 
symmetric with respect to interchanging a and /3-sites. 
Since we precluded the possibility of supercritically en- 
hanced surface bonds, a spontaneous breakdown of this 
symmetry is ruled out for T > T^. Therefore the solu- 
tions that minimize the free energy in this temperature 
regime fulfill m„_Q, — m„^^ = m„, and G acts in a 2D sub- 
set {{x,x,y,yY'\x & \y\ < 1} of R''. The picture of 
the invariant manifolds looks similar to Fig. 5. The mag- 
netization profile at the bulk critical point K = Kc{h) 
(for /i = 1) is depicted in Fig. 7. As in the case of the 
(100) orientation, the layer magnetization densities oscil- 
late about the bulk value mdis due to the antiferromag- 
netic coupling between neighboring planes. However, the 
decay length remains on the order of the lattice constant 
even at T = Tc, (cf. Sec. Ill A). The numerical profile 
decays exponentially on a length scale that agrees well 
with the value (36b). 
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FIG. 6. (a) Magnetization profile nin for the (100) surface 
above Tc, as obtained from the intersection S of Fig. 5. Odd 
and even planes n belong to sublattice a (0) and /? (x), re- 
spectively. Note that mn oscillates about the bulk value mdis 
(dashed line), (b) OP profile, defined by (32). For conve- 
nience, the discrete points have been joined by lines. 



3. (100) surface, T < Tc 

On crossing the critical line in the MF phase diagram 
(Fig. 2), the map F undergoes a period-doubling bifur- 
cation (Sec. Ill A), and the picture of the invariant mani- 
folds changes qualitatively. The fixed point Vdis becomes 
elliptic and looses any inset and outset. The inset (out- 
set) of any one of the hyperbolic fixed points v^f^ cannot 
intersect itself but will generically intersect the outset 
(inset) of the same fixed point, as well as the outsets 
(insets) of all other fixed points, at an infinite number 
of so-called homoclinic and heteroclinic points. The in- 
variant manifolds oscillate wildly in the vicinity of the 
fixed points (cf. Fig. 8), giving rise to the phenomenon 
of "chaotic entanglement".^'^ As a consequence, infinitely 
many solutions of the lattice MF equations for the semi- 
infinite system exist, and one has to resort to the original 
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FIG. 7. (a) Magnetization profile at T = for tlie (110) 
surface. The dashed line represents the bulk value rridis. (b) 
Exponential decay of the profile. The solid curve is a straight 
line with slope [Eq. (36b)]. 



variational principle (cf. Sec. II) in order to decide which 
solution corresponds to the true equilibrium profile. The 
stable and unstable manifolds Wg and Wu of v\^^ and 
Dqj.^ are determined uniquely given only one of them, 
say Ws{y\^^ . In fact, from the symmetry property (19) 
one concludes that 

Wu{vl^) =R{Ws{vl^)). 

Typical minimum-free-energy profiles are shown in Fig. 
9a. The corresponding trajectories converge to the fixed 
point Vi describing the bulk phase 0^ > 0. If one imposes 
the bulk boundary condition </>„ ^ — 0b < (n — > cxd), 
so that the solutions converge to V2, one obtains OP pro- 



files exhibiting an antiphase boundary, i. e., an interface 
between the two phases ±(/)b (Fig. 9b). These solutions 
always yield a higher free energy than the profiles of Fig. 
9a. For the chosen parameter values, the free energy of 
the profiles is found to increase as the position of the 
interface moves into the bulk, so that the leftmost profile 
represents the equilibrium solution for this type of bound- 
ary conditions. By analogy with wetting phenomena,^" 
we may say that the surface is "nonwet", i. e., the inter- 
face has a finite (microscopic) distance from the surface. 
If hi is increased, so that the "effective" ordering sur- 
face field gi favoring the bulk phase 0;, > becomes 
strong enough, the free energy of the profiles eventually 
decreases and the surface is "wet" . (E. g., in the situation 
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FIG. 8. (a) Part of the stable manifold of the fixed point 
v^j.^. The inset shows the first few loops in the vicinity of 
''^ord- (b) Magnification of the region marked in the inset of 
(a). The dashed-dotted line is the boundary condition for 
hi = 0.39, which intersects the stable manifold at an infinite 
number of points Si, S2, .... The latter accumulate at the 
point Soo on the stable manifold of uJrd (dashed line). 
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FIG. 9. (a) OP profiles below Tc for the (100) surface. The 
corresponding trajectories converge to the fixed point v^^^ (cf. 
Fig. 8b). The middle profile {hi = 0.39) belongs to the inter- 
section Soo- (b) OP profiles describing an antiphase bound- 
ary. The trajectories converge to v^^.^. The leftmost profile 
belongs to the intersection 5*1, the next one to 5*2, etc. (Fig. 
8b). An infinite number of such solutions exists, the first ten 
of which are shown here. 



boundary condition for the OP profile (/^{z), 

where z = corresponds to the surface plane (n — 1) and 
the dot denotes differentiation with respect to z. Such a 
boundary condition is familiar from the phenomenologi- 
cal theory of surface critical phenomena."^ The parame- 
ter C > is the coefficient of the gradient term of the 
Ginzburg-Landau functional. The extrapolation length 
X should be positive owing to the absence of enhanced 
surface bonds. Thus the persistence of surface order for 
T > Tc originates solely from the "effective" ordering 
surface field gi ^ 0. We will determine the dependence 
of gi — gi{K, h, hi) on the reduced coupling constant K 
and the fields h and hi. 

No such ordering surface field emerges in the case of 
the symmetry-preserving (110) surface. However, the 
Ginzburg-Landau functional depends on an additional 
spatially varying nonordering density to which the sur- 
face field hi couples. We defer the derivation of the cor- 
responding continuum model to a subsequent paper, 
where nonordering densities and the construction of suit- 
able multicomponent Ginzburg-Landau theories will be 
treated from a more general perspective. 

Upon approaching the bulk critical point in the pres- 
ence of an ordering surface field gi ^ 0, the semi-infinite 
system undergoes the normal transition, which exhibits 
critical singularities that are distinct from those of the 
ordinary transition. The latter occurs for gi — Q and 
subcritical surface enhancement (A > 0).^'^ In order to 
confirm that the continuum theory derived in Sec. IV A 
correctly describes the asymptotic critical behavior of the 
lattice model, we shall draw a detailed numerical com- 
parison with the solutions of the lattice MF equations 
and clearly identify the singular behavior of the normal 
transition for generic values of h and hi (Sec. IV B). By 
tuning h and /ii, one can achieve that gi = Q aX T = Tc- 
Then the singularities of the ordinary transition are re- 
covered, although the OP at the surface is nonzero for 
T>Tc because oi gi ^ Q for T ^ Tc (Sec. IV C). 



depicted in Fig. 9b this happens if one chooses hi — 0.4.) 
In Ref. 21 the wetting phase diagram in the space of ther- 
modynamic parameters K , h, and hi is calculated using 
the continuum model derived in the next section. 



IV. GINZBURG-LANDAU THEORY 

In this section, the aim is to derive and critically exam- 
ine a Ginzburg-Landau model for the semi-infinite alloy 
with a (100) surface. In particular, we want to show that 
the loss of the a ^ /3 sublattice symmetry (cf. Sec. I) 
leads, in a continuum description, to a symmetry-breaking 



A. Derivation of continuum model 

One has to be careful to define the local OP (pn in 
such a way that the space-inversion symmetry of the lat- 
tice MF equations (8) survives the continuum limit. If 
the magnetization profile m„, n = 1,2,..., is a solu- 
tion of (8), so is the profile m_„ = to„+2 obtained by 
reflection at the surface plane n = 1. If the definition 
of (j)n respects this symmetry, the differential equation 
for the continuum profile (j){z) should be invariant under 
4>{z) (j){—z). However, the definition (32) distinguishes 
one direction along the [100] axis and violates space- 
inversion symmetry explicitly. As a consequence, first- 
order derivatives of (l){z) appear in the Ginzburg-Landau 
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equations, corresponding to linear derivative terms in the 
free-energy functional.^ Such terms render the functional 
unbounded from below and thus preclude it from serv- 
ing as a Landau- Ginzburg- Wilson Hamiltonian of a field 
theory. 

To avoid these difficulties, we adopt an alternate defi- 
nition of the local OP treating the two neighboring layers 
of lattice plane n on an equal footing: 



-i-lY 



1/ 

-{mn-i+m„+i) 



rrir, 



(37) 



This definition complies with space-inversion symmetry 
and coincides with (32) in the bulk of the system. 

Substituting TO„+i-|-m„_i from (8) into (37) we obtain 



= (-1)"+V(mn) , 



(38) 



where the function 

fix) = 



h 

4 



X 

2 



4:K 



■ tanh 



is strictly monotonous and thus invertible. We denote 
the inverse of by , so that from (38) 



m„ =M((-1)"+Vn) 



(38') 



The continuum limit of the lattice MF equations (8) leads 
to the Ginzburg-Landau equation (see Appendix B) 



Co^ = 40 + 2Al(-0) - 27W(0) , 



(39) 



where C = M'{0) [Eq. (B4a)]. The MF equation at the 
surface (9) entails the boundary condition (see Appendix 
C) 



C^(0) = -hi- Mi-(p{0)). 



(40) 



For a spatially homogeneous system, (39) is identical to 
the equation for the OP following from the bulk MF equa- 
tions (A2'). In fact, (A2') can be written as 

(f) = tp{m + <!)), -(p = <^(m - (j)) . 

Operating with A4 on both sides of the above equations 
and eliminating m, one arrives at (39) with ^ = 0. 

We denote the bulk Landau free-energy density by 
fb{ma,m0) [Eq. (Al)]. Since 

dafbirria, mp) = 2ip{ma) - {rua - m^) , 
dpfbirria, mp) = 2ip{m0) + (m„ - m^) , 

where 9^ = fi = a,P, Eqs. (39) and (40) can be 

rewritten as 



C4> = v' {<!>), 

cm = fs{<p) , 



(39') 
(40') 



where 



Vi<f>) = -[M{cP)-M{-(t>)-2ct>f 

+ h{M{<j>)M{-<i^)) - h{M{QlMm , (41) 
fsi<f>) = -[hi + M{-<t>M+\MMi-ct>),M{-<t>))- 



--MM{0),Mm. 



(42) 



Thus the Ginzburg-Landau equations (39) and (40) fol- 
low from the variation of the free-energy functional 



m = / 

Jo 



dz 



+fs{m)- (43) 



Expansion of V{(f)) yields the usual (p^-foim 

n0) = |^' + f^'+o(/), 



(44) 



where ^4 = 4(1 - Mi) and B = -AMz [Eqs. (B4)]. With 
the aid of (A7), the leading temperature dependence of 
the Landau coefficients is found to be 



A = Ait + 0{f) , 
B = Bo + Oit) , 
C = l + Oit) . 



(44a) 
(44b) 
(44c) 



where 



Al = Ai{h) = 2{1 - Kcih)mc{h)[h - 2toc(/i)]} , 
Bo = Bo{h) = ^i^c(/i)' . 

Likewise, the Landau expansion of fs{<p) reads 

C 



where 



9i{K, h,hi) = hi+ mdis{K, h) , 
A = l. 



(45) 



(45a) 
(45b) 



Thus we obtain a positive extrapolation length A and an 
ordering surface field gi ^ 0, as anticipated above. 



B. Comparison with lattice MF theory: generic 
(nonideal) stoichiometry 

If h and hi take generic values (nonideal bulk and 
surface stoichiometry), gi = gi{K,h,hi) is nonzero at 
T = Tc and gives rise to an OP profile decaying accord- 
ing to a power law. For t = 0, (39') becomes 



(46) 



The neglected higher powers of (p do not affect the asymp- 
totic behavior of (piz) as z ^ 00. The solution of (46) 
satisfying (f){z) — > as 2: — > 00 reads 



11 





1.74 




1 1 1 




1 79 




/i = 1 , K = Kc{h) - 




1.7 




hi = 0.4 


p 


1.68 






-1- 

1 


1.66 
1.64 


o 
-o 

o 
- o 

% 






1.62 
1.6 








1 1 1 1 



20 40 60 80 100 

lattice plane n 



FIG. 10. Scaling of the critical OP profile [defined by 
(37)]. The expected power law decay has been factored out. 
The parameter no was fitted optimally. The dashed line is 
the prediction (46b) of the continuum theory. 

Hz) = , (46a) 

Z + Zq 

with the amplitude 

The signs refer to positive and negative gi, respectively. 
The integration constant zq follows by inserting (46a) 
into the boundary condition (40). An algebraic decay 
4>{z) ~ 2;"^* of the OP profile at T = Tc, where is 
the scaling dimension of the OP, is characteristic of the 
normal transition. Generally, a;^ can be expressed by 
bulk critical exponents as = P/v. Within MF theory, 
/3 = ^ = 1/2, so that Xcj, — 1. 

In order to check whether the lattice MF profiles ex- 
hibit the above power law decay with the predicted value 
(46b) for the amplitude, we used the "nonlinear-mapping 
representation" of the lattice MF equations (Sec. Ill) to 
determine the profiles numerically up to ~ 1000 layers 
from the surface. Fig. 10 demonstrates that the power 
law decay is well reproduced. The amplitude extracted 
from the numerical fits is in perfect agreement with (46b) 
for a wide range of the bulk field h (Fig. 11). 

As a second test of the validity of our continuum model 
we investigate the temperature singularity of the surface 
OP. We multiply the Ginzburg-Landau equation (39') by 
(j) and perform the integral from z = to z = oo on both 
sides using the boundary conditions (40') and 4>{z) 
as z — > cx). This gives 

~^f'si'l^sr + V{c^s) = V^, (47) 

where (ps = ^(0) and Voo = V{<j)b)- One has 




magnetic field h 

FIG. 11. Amplitude of the critical OP profile as a function 
of the bulk field h. The data obtained from the solution of 
the lattice MF equations (□) agree well with the prediction 
(46b) of the continuum model (full line). 



^^10, t>0, 

[-At' + o{t'), t<0. 

Expanding the coefficients on both sides of the surface 
equation of state (47) in powers of t, one recognizes that 
(j)s exhibits a discontinuity in the second temperature 
derivative due to the non-analyticity of Voo at t = 0, 

0.-4°)+0Wt + 0i'it2 + ... , (48) 

with (pl ^ (pl '_. This result complies with the general 
form of the leading singularity of the surface OP 

at the normal and extraordinary transitions.^^ 

Fig. 12 shows a comparison with the solutions of the 
lattice MF equations. While the singularity in the sec- 
ond temperature derivative is too weak to be detected 
without considerable numerical effort, the data are fully 
consistent with the continuity of the first temperature 
derivative at i = and agree well with the predictions of 
the continuum theory. The continuity of the first temper- 
ature derivative of (ps is incompatible with the ordinary 
transition and confirms again that the asymptotic criti- 
cal behavior falls into the universality class of the normal 
transition. 



C. Comparison with lattice MF theory: Vanishing 
ordering surface field at T = Tc 

According to (45a) , gi can be made to vanish at T = Tc 
by choosing hi such that 

hi = -mc{h) . (49) 
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FIG. 12. Temperature dependence of the surface OP <l)s for 
generic values of h, hi , in the immediate vicinity of the critical 
temperature. The results of the lattice MF theory (□) agree 
well with the prediction of the continuum model (solid line). 



If h = hi = (ideal bulk and surface stoichiometry), gi = 
for all temperatures, and the system clearly displays 
ordinary surface critical behavior, li h ^ and (49) is 
fulfilled, gi varies linearly with t as t ^ 0. In particular, 
the surface OP <f>s is nonzero for T > and vanishes 
only in the limit t — > 0±. One may wonder whether 
such a behavior is consistent with the ordinary transition 
where one usually expects that 0s = for T > Tc- 

To derive the leading temperature singularity of (pg 
consider again Eq. (47). Since (/)s = if t = 0, 0s is 
found to behave asymptotically as 



±t + less singular terms , 



(50) 



with 



7^ i^i^l. The ordering surface field (45a) can 



be expanded as 

gi{K,h,hi)=g^^\h,h,)+g^^\h)t + 0{t^) , (51) 
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FIG. 13. Temperature dependence of the surface OP in the 
case where the effective ordering surface field gi vanishes at 
T — Tc- The discontinuity in the first temperature derivative 
predicted by the continuum model (full line) is well confirmed 
by the numerical solutions of the lattice MF equations (□). 



of the lattice MF equations (Fig. 13). As will be shown 
below, such a behavior is precisely what one expects if 
the leading asymptotic behavior belongs to the univer- 
sality class of the ordinary transition. The variation of 
the effective ordering surface field gi with temperature 
[cf. Eq. (51)] explains the onset of surface order for t > 
(see below). 

Let us elucidate the above behavior of the surface OP 
05 by resort to a scaling argument. The basic assump- 
tion is the existence of a scaling field gi = gi{K,h,hi) 
associated with 0^ and depending analytically on K, h, 
and hi. Of course, gi will in general differ from the MF 
expression (45a). By analogy with (51) we write 



gi{K,h,hi 



g^i\h,hi) 



g['\h,hi 



(54) 



We suppose that gf \h, hi) vanishes if h and hi are cho- 
sen appropriately. This should always be possible since 
where g[^''{h,hi) = hi + mdh) = owing to (49) and 0s is positive for large positive hi, and becomes negative 

if one lets hi — > — oo, assuming h and t to be fixed. Thus 
for given h and t — 0, 0^ must vanish for a special value 
of hi, in which case one clearly has ^J"' = 0. 

The singular part of the surface free energy /s''"^ 
should take the standard scaling form'^ 



g[^\h) = TOd2(/i) [Eq. (A7)]. Insertion into (47) yields, 
to leading order in t. 



So 



t > 0, 
t < 0. 



Therefore one concludes that^^ 

- .(1) 



9i 



Ai 



(52) 



(53) 



fr%t,~gi)^Mi'-'^''\t\^^-'^'^g±[ 



Mg,M-^^~gi\t\ 



-Ai 



The discontinuity of 0^ in the first temperature deriva- 
tive differs strikingly from the singularity in the second 
derivative for generic values of h, hi (Sec. IV B). This re- 
sult is in excellent agreement with the numerical solutions 



where Ai ~ 0.48. All non-universality is embodied in 
the metric factors Mt and Mg^ associated with the two 
relevant scaling fields (at the ordinary transition) t and 
gi , while the critical exponents and the scaling functions 
5±(C) are universal. The singular part of 0s follows by 
taking the derivative with respect to gi, i. e.. 



^A'CMg,\tf'Y±{Mgi\t 



(55) 
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where Y± = g^, M = Mg^Mt^\ and 01 = {d - l)v - 

Ai. The regular contribution cji^s^^^ describes, to leading 
order, the linear response of 



(56) 



The importance of such regular terms for the correct 
identification of surface critical exponents and scaling 
functions has been emphasized in Rcf. 26. The scaling 
functions 1±(C) are analytic at ( = 0, 



(57a) 
(57b) 



Y+{C)=Yi'\ + 0{e) , 

y_(C) = Fi°^ + yi'V + o(C') • 

Using (54) with g[°^ = 0, Eqs. (55)-(57) yield 

<t>s = M,^^M-g,\tf^Y^ (Mg['k\t\-^^) +ag['k + 0{t^) , 

where terms of order f^'jij^^^, p = 2,3,..., have been 
omitted in the argument of Y±{(). Thus the leading be- 
havior as t ^ becomes 



ag[^h + M^' Mg,Mg^^^Y!^^h^- 
M^'M~g,Y^'^\t\l^^ 



■71,1 



t>Q, 
t<0. 



(58) 



where wc used the scaling relation —71.1 = /3i — Ai. Since 
7i,i < at the ordinary transition,"^ 0^ varies linearly 
with f as t — > 0+, but vanishes with the characteristic 
exponent /3i ~ 0.8 of the ordinary transition as t ^ 
0— . In MF theory, /3i = 1 and the power singularity for 
t ^ 0— degenerates into an integer power, see Eqs. (50) 
and (53). That the asymptotic behavior of the ordinary 
transition can be obtained by tuning h and hi has also 
been demonstrated by transfer matrix calculations in two 
dimensions,^'' which supplement the results of Ref. 8. 



V. SUMMARY 

We have studied the surface critical behavior of bcc 
binary alloys undergoing a continuous B2-A2 order- 
disorder transition in the bulk. Clear evidence has been 
found that symmetry- breaking surfaces, such as the (100) 
surface, generically display the critical behavior of the 
normal transition, which belongs to the same universality 
class as the extraordinary transition. Wc have analyzed 
the lattice MF equations using the "nonlinear-mapping" 
representation^^ and achieved a mapping onto a contin- 
uum (Ginzburg-Landau) model. The latter assumes the 
form of the standard one-component (^^-model for semi- 
infinite systems. Its crucial feature is the emergence of an 
"effective" ordering surface field gi 0, which depends 
on temperature and the other parameters of the lattice 
model and is not present on a microscopic level. By a 



detailed comparison with the solutions of the lattice MF 
equations the continuum model has been shown to ac- 
curately describe the asymptotic behavior of the lattice 
model. 

In the case of the symmetry-preserving (110) surface 
the appearance of an ordering surface field is ruled out by 
symmetry. Analysis of the lattice MF equations reveals 
the existence of an additional length scale different from 
the OP correlation length, which describes the decay of 
the nonzero magnetization profile above T^. However, 
this length stays microscopic even at T = Tc and does 
not infiuence the leading singular behavior, which is char- 
acteristic of the ordinary transition. The construction of 
a suitable continuum model in this case is deferred to 
a subsequent paper, where particular emphasis will be 
laid on nonordering densities and the derivation of the 
corresponding multicomponent Ginzburg-Landau theo- 
ries from microscopic models. 
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APPENDIX A: BULK MF EQUATIONS 

For a spatially homogeneous system with sublattice 
magnetization densities rria and m^, the variational free 
energy (6) yields the Landau free-energy density 

„ . X Fmfa/N h 

h[ma,m/3) = — — = niamfs - -{rUa + m/s) 



+ 



1 

2K 



dx tanh" 



1 

Jo 



dx tanh ^ a; J , (AI) 



where is the number of lattice sites. The MF equations 
dafb = 0, d0fb = O, (A2) 



where 9^ = = a, read 

tanh~^ rria = K{h — 2to^) , 
tanh""^ m/3 = K{h — 2ma) ■ 

For the following it is useful to define 

Ua = U{mct) , Uff = u{mi3) , 

where 

1 1 



u{x) 



2Kl-x^' 



(A2a) 
(A2b) 



(A3) 



(A4) 
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A solution {ma,mp) of (A2) minimizes fb and is thus 
thermodynamically stable if the matrix 



dafb dadpfb 



dadpfb djfb 
is positive definite, i. e., 

UaU/3 > 1 . 




(A5) 



If one writes nia = m + cj) and rn/3 = m — </> as in (21), 
the MF equations (A2) take the form 

tanh"^ (m + (!>) = K{h - 2m) + 2K(t> , (A2a') 
tanh"^ (m - (t>) = K{h - 2m) - 2K(l) . (A2b') 

The disordered state (m = mais, = 0) satisfies 

tanh""^ mdis = K{h — 2mdis) • (A6) 

This equation has a unique solution Wdis = iTidis{K,h), 
which may be expanded in powers of the reduced tem- 
perature t = (if c — K)/K as 

mdis = mc + m^2* + ^ if) > (^7) 
where rric = mc{h) is the magnetization at T = T^, and 



(A7a) 



The disordered state is thermodynamically stable only if 
w(mdis) > 1 [cf. Eq. (A5)]. The phase transition occurs 
when u(mdis) = 1 [cf. Eq. (20)]. Two new minima of 
fb describing the ordered phases (m = mord,'/' = i^fc) 
emerge if w(TOdis) < 1. The asymptotic behavior follow- 
ing from (A2') is found to be 



Word = TOc + m'^^^t + O (f) , 



(A8) 



where 



= ml^J^ih) = m2{h) + K,(h)m,{h)Uhf , (A8a) 



V3 



2K^{h) 



V'l - Kc{h)mc{h)[h - 2mc{h)] . 

(A8b) 



In the continuum limit, one replaces 4>n by a smooth 
profile (t>{z) defined for all z > 0, with the original layers 
located at z„ = n — 1. Assuming that the OP varies 
slowly on the scale of the layer spacing, we approximate 



M{{-ir<t>n-l)+M{{-ir<f>n+l) 

~ 2M{i-ir<f>izn)) + i-irM'mizr,) , 



(B2) 



where the dot denotes differentiation with respect to z 

and terms of order 4>(j) and (j)^ have been discarded. Sub- 
stitution of (B2) into (Bl) leads to the Ginzburg-Landau 
equation (39). 

Since (^(mdis) = 0, i. e., A^(0) = mdis, A^(0) may be 

expanded as 

M{ct>) = mdis + Mi0 + M2</.2 + M3</.3 + O (</>4) , (B3) 
where 



Ml =7W'(0) = 2[l + u(mdis)]-S 
M2 = -Km^iMmAisfMl , 

Ms = ^{sMi [2ifmdisw(TOdis)^]' 



(B4a) 
(B4b) 



- {2Kf (1 + 3mD u{maisf}Mt . (B4c) 



APPENDIX C: BOUNDARY CONDITION 

By analogy with (Bl), the MF equation at the surface 
(9) can be written as 



The continuum approximation (B2) now reads 



X(-02) ^Mi- </.(0)) - A^'(O) 



(Cl) 



(C2) 



Inserting (C2) into (Cl) and requiring that the Ginzburg- 
Landau equation (39) be also valid at 2 = 0, we obtain 

hi + M{- (l>{0)) + M'{0)^{0) = ^M'{0)^{0) . 

The above equation reduces to the boundary condition 
(40), if second-order derivatives are neglected. 



APPENDIX B: GINZBURG-LANDAU EQUATION 

Using (38) and (38'), we may rewrite the MF equations 
(8) as 

A^((-l)"</.„_i) +M{{-ircf>n+i) 
= 2A^((-l)"+Vn)+4(-l)".^„. 



^ H. Dosch, Critical Phenomena at Surfaces and Interfaces, 
Springer Tracts in Modern Physics Vol. 126 (Springer, 
Berlin, 1992). 

^D. M. Kroll and G. Gomppcr, Phys. Rev. B 36, 7078 
(1987); G. Gompper and D. M. Kroll, Phys. Rev. B 38, 
459 (1988). 



15 



^ K. Binder, in Phase Transitions and Critical Phenomena, 
edited by C. Domb and J. L. Lcbowitz (Academic, London, 
1983), Vol. 8, p. 1. 

* H. W. Diehl, in Phase Transitions and Critical Phenomena, 
edited by C. Domb and J. L. Lebowitz (Academic, London, 
1986), Vol. 10, p. 75. 

^ H. W. Diehl, in Proceedings of the Third International 
Conference "Renormalization Group 96", Dubna 1996, 
to be published, and preprint cond-mat/9610143. 

^ F. Schmid, Z. Phys. B 91, 77 (1993). 
The crucial feature of both transition is that the bulk or- 
ders (at T = Tc) in the presence of an already ordered 
surface. The extraordinary transition exists only in bulk 
dimension d > 2 (in the Ising case), i. e., if the dimension 
of the surface is greater than the lower critical dimension. It 
can be shown that the asymptotic surface critical behav- 
ior is the same, irrespective of whether the surface order 
occurs spontaneously due to supercritically enhanced sur- 
face interactions (extraordinary transition), or is induced 
by an ordering surface field (normal transition), see T. W. 
Burkhardt and H. W. Diehl, Phys. Rev. B 50, 3894 (1994). 

® A. Drewitz, R. Leidl, T. W. Burkhardt, and H. W. Diehl 
Phys. Rev. Lett. 78, 1090 (1997). 

^ S. Krimmol, W. Donnor, B. Nickel, H. Dosch, C. Sutter, 
and G. Griibel, Phys. Rev. Lett. 78, 3880 (1997). 
The width of the asymptotic region can be estimated by 
requiring that the scaling variable be on the order 

of one, where t = (T - Tc)/Tc and Ai ~ 0.48. If we esti- 
mate the magnitude of gi using the MF expression derived 
below [Eq. (45a)], we conclude that normal critical behav- 
ior should be observed if \t\ < 0.06. However, the closest 
approach to Tc in the simulations of Ref. 6 was \t\ ~ 0.14. 
" R. Pandit and M. Wortis, Phys. Rev. B 25, 3226 (1982). 
F. Ducastelle, Order and Phase Stability in Alloys, (North 
Holland, Amsterdam, 1991). 

D. K. Arrowsmith and C. M. Place, An introduction to dy- 
namical systems (Cambridge University Press, Cambridge 
UK, 1990). 

"'^ The reentrant behavior for \h\ > 2, i. e. the possibility 
of passing from the ordered to the disordered phcise upon 
lowering the temperature at fixed h, is not confirmed by 
Monte Carlo simulations and cluster-variation calculations, 
see Ref. 15. Therefore, we will always assume that \h\ < 2 
in the following, so that T < Tc {T > Tc) corresponds to 
the ordered (disordered) bulk phase. 
B. Diinweg and K. Binder, Phys. Rev. B 36, 6935 (1987). 
As we shall see in Sec. IV, a slightly different definition of 
the local OP is favorable in order to perform the continuum 
limit of the lattice MF equations. However, this alternate 
definition changes none of the subsequent conclusions. 
V. Privman, P. C. Hohenberg, and A. Aharony, in Phase 
Transitions and Critical Phenomena, edited by C. Domb 
and J. L. Lebowitz (Academic, London, 1991), Vol. 14, p. 
1. 

1* The lengths and i [Eqs. (31a) and (34a)] differ by a 
factor of \/2 (to leading order as t — > 0) since we implicitly 
used the layer spacing as a unit length. Therefore the bulk 
correlation length is y/2 times larger when measured in 
units of the (100) rather than the (110) layer spacing. 
Since the fixed point corresponding to the thermodynam- 



ically stable bulk phase is unstable under the nonlinear 
mapping, any trajectory computed numerically will even- 
tually flow away from the fixed point, but the error in- 
volved when calculating, e. g., order parameter profiles can 
be made arbitrarily small using standard techniques. 
S. Dietrich, in Phase Transitions and Critical Phenomena, 
edited by C. Domb and J. L. Lebowitz (Academic, London, 
1988), Vol. 12, p. 1. 

R. Leidl, A. Drewitz, and H. W. Diehl, preprint cond- 

mat/9704125. 

R. Leidl, to be published. 

For the ordinary-normal crossover in the presence of a small 
gi and a discussion of the anomalous short-distance behav- 
ior of the OP profile see Ref. 5, pp. 17-19, and A. Ciach 
and U. Ritschel, Nucl. Phys. B 489, 653 (1997). 
Beyond MF theory, the amplitude ratio \ of the 

singularity of 0s is known to be identical to the uni- 
versal value B+ /B- , where B+ and B- are the amplitudes 
of the " singularity of the bulk free energy (Ref 7). 
However, in MF theory one cannot separate the singular 
contributions in the expansion of (bs from the regular back- 
ground terms. Therefore, the ratio (f>l \./4'l _ is a nonuni- 
versal quantity in MF theory. 

Without loss of generality we may assume that h > 0, so 
that gi is positive for T < Tc. Then the second solution of 
(47) for t < 0, which would lead to cp'^l = g[^^ + Ai/y/B^, 
yields a higher surface free energy and may thus be dis- 
carded. 

A. J. Bray and M. A. Moore, J. Phys. A 10, 1927 (1977) 
A. Drewitz, unpublished. 



16 



